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We study the scattering of graphene quasiparticles by topological defects, represented by holes, 
pentagons and heptagons. For the case of holes, we obtain the phase shift and found that at low 
concentration they appear to be irrelevant for the electron transport, giving a negligible contri- 
bution to the resistivity. Whenever pentagons are introduced into the lattice and the fermionic 
current is constrained to move near one of them we realize that such a current is scattered with 
an angle that depends only on the number of pentagons and on the side the current taken. Such 
a deviation may be determined by means of a Young-type experiment, through the interference 
pattern between the two current branches scattered by a pentagon. In the case of a heptagon such 
a current is also scattered but it diverges from the defect, preventing a interference between two 
beams of current for the same heptagon. 
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1. Introduction and Motivation 

Condensed matter physics is a branch of physics which studies systems of many particles in 
the condensed, i.e. solid or liquids states. The Schrodinger equation is the start point of current 
condensed matter theory, however, in practice an generic interacting many-body system is an ex- 
tremely complex system and it may not always be helpful for obtaining physical properties of an 
interacting many body system. The collective excitations of a many body system are collective 
motions of the atoms and can be viewed as particles, however, the properties of those particles can 
be very different from the properties of the particles that form the many body system. 

An example of this collective behavior is the graphene a flat monolayer of carbon atoms tightly 
packed into a two-dimensional (2D) honeycomb lattice, consisting of two interpenetrating triangu- 
lar sublattices [|J]. It is the first example of a truly atomic two-dimensional (2D) crystalline system 
and it is the basic building blocks for graphitic materials such as fullerenes (graphene balled into a 
sphere) or carbon nanotubes (graphene rolled-up in cylinders) [0, |J]. Experimental techniques 
provide high-quality graphene crystallites, up to 100/im in size, which is sufficient for most re- 
search purpose, including the ones considered here. In this "perfect" layer the charge carries can 
travel thousands of interatomic distances without scattering. 

Graphene provides an bridge between condensed matter physics and quantum electrodynamics 
because the colletive excitations are described by Dirac equation for massless particles in (2+l)D, 
it is a zero-gap semiconductor, in which the low energy spectrum is correctly described by the 
(2 + 1 )D Dirac-like equation for a massless particle []|] 

ft— |¥)=v F ff-/f|¥), (1-1) 

where vp is the Fermi velocity, which plays the role of the speed of light (vp ~ c/300), a = (o x , o y ) 
are the 2D Pauli matrices, p = —ihV is the linear momentum operator and is a two-component 
spinor. Therefore, the quasiparticles can be viewed as electrons that have lost their masses or as 
(massless) neutrinos that acquired the electronic charge. Such a spectrum makes graphene a ma- 
terial with unique electronic properties. Its description by means of the Dirac equation is a direct 
consequence of graphene's crystal symmetry. Its honeycomb lattice is made up of two equivalent 
triangular carbon sublattices A and B, whereas its cosine-like energy bands associated with the sub- 
lattices intersect at zero energy (E = 0) near the edges of the Brillouin zone, giving rise to a conical 
section spectrum at low energies^], say, \E\ < 1 eV. In this honeycomb lattice, the two-component 
spinor |*F) is referred to as pseudospin since it is an index indicating two interpenetrating triangular 
sublattices AeB, which is similar to spin index (up and down) in quantum electrodynamics (QED). 
It is common to regard the sublattice degree of freedom as a pseudospin, with the A sublattice be- 
ing the "up", |+) and B sublattice being the "down", — }. Since vf <C c, it is a slow relativistic 
system or a strong coupling version of QED since the graphene's dimensionless coupling constant, 
e 2 /fivp m 1 much higher than its QED analogue, the fine structure constant e 2 /Tic « 1/137. All 
these properties make the graphene a very interesting system, which provides a way to probe QED 
phenomena, for instance, by measuring its electronic properties. Several proposals for testing some 
predicted, but not yet observed phenomena in QED, including the Klein paradox [j5], |6|], vacuum 
polarization [f7|] and atomic collapse are some topics under investigation in graphene. 
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Here, we would like to study the behavior of graphene quasiparticles in the presence of de- 
fects in the crystalline structure of the material. We shall consider three types of defects: holes, 
pentagons and heptagons [§J. All these defects can be incorporated by removing or inserting a 
few carbon atoms. The presence of defects like pentagons (heptagons) induces positive (negative) 
curvature in the material. At some extent the charge carriers motion in the presence of pentag- 
onal (heptagonal) defects is identical to fermions moving in a (2+ 1)D gravitational space-time 
generated by positive ("negative") point-like masses. Then we may employ results from gravity 
to analyze some effects concerning charge carriers in the presence of such defects. Understanding 
how these defects modify the transport properties of graphene is crucial to achieve future electronic 
devices using carbon-made materials. 



2. Scattering of graphene quasiparticles by holes 

In the continuum model for the graphene, it is assumed that there is a hole of radius r$ cut 
from the system center, located at the origin. Thus the motion of the quasiparticles is performed 
on a flat 2D support given by a non-simply connected manifold, which can be viewed as defect in 
the material. This model allows for the investigation of scattering effects as a function of the hole 
radius and could shed some light on the high charge carrier mobilities, fact observed in graphene 



[10] 



To determine the phase-shift of the scattered wave function as well as the scattering cross 
section one has to solve the two-dimensional Dirac equation ( |1.1| ) which, for the case of massless 
particles, can be writing in a covariant form, ihy^d^\j/(x) = 0, where the covariant derivative is 
dfi = [(l/vp)d/dt , d/dx, d/dy], the y-matrices are y° = a 3 , y 1 = io 2 and y 3 = — ia , obeying 
yfiyv _ rjiiv _j £ [ivay a ^ is the Minkowski tensor metric, diag(r]' IV ) = (+1,-1,-1) and e^ va 
is the 3-dimensional Levi-Civita symbol (e 012 = +1). (The word covariant must be used carefully 
because vp is not invariant, being only a parameter and the term covariant is used to refer only the 
form the equation is written.) We may expand the solutions of the free massless Dirac equation in 
plane waves, once rotational invariance allows to separate the 6 variable, so that the diagonalized 
angular momentum J? = —iftjg + §<7 3 , yielding partial waves with angular momentum (n + \)fi, 
takes the form, yr(r,t) = e i(n+ ^-^ e u n {r)e- iE 'l h . 

The components of the radial spinor u n (r), given by /„ (r) and g„ (r), satisfy the Bessel equation 
of order n and n + 1. The solution for the radial spinor outside of the hole (i.e., for r > ro) is given 
by: 

(kr)+B 2n N n (kr) \ 
U " [l) 1 ! " \B 3n J n+l (kr)+B 4n N n+l (kr) )' 

where J n and N n are the Bessel functions of first and second kinds (Neumann function), respec- 
tively, n = , ±1 , ±2 , ... is the angular-momentum number and k = ^ > (we are considering 
only solutions with E > that describe the electronic dynamics) and Bj n (j = 1 , 2 , 3 , 4) are con- 
stants. These constants must be determined by the appropriate boundary conditions specified to 
completely define the problem. From the physical point of view, the correct boundary condition 
is determined by the requirement of vanishing the net energy flux into the hole, which is a region 
absent of lattice degrees of freedom. Consequently, the fields must arrange themselves in such a 
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way that the energy flux from the incoming modes (asymptotically behaving like e~' kr ) exactly 
cancel that from the outgoing waves. Imposing Neumann boundary condition (NBC) on the wave- 
functions, ^ 



= 0, we obtain: 

r=r 



u ( r )=( B ln [j n (kr)-tan{t n (kr ))N n (kr)) \ 
Un{n \B 3n [j n+l (kr)-t a n(t n+l (kr ))N n+l (kr)] J 

where tan[t n (kro)} = jfr^p^j • The terms proportional to Bessel (Neumann) functions describe the 
incident (scattered) waves. Comparing their asymptotic behavior with those for plane waves, giving 
us the phase-shift 8 n of the n-th partial wave that completely determines the fermionic scattering: 

8*=t n (kr ), S* = t n+1 (kr Q ). (2.3) 

Now, if we consider a small concentration of point-like defects with concentration n^ef angle- 
depedent scattering cross section, o(6), their contribution to resistivity, p, may be estimated as 
P — «def -jfr - (kro) 6 [§]. This means that the scattering induced by small holes (with radius around a 
few angstroms; some lattice spacings) at low concentration are irrelevant for the electronic transport 
in graphene, giving a negligible contribution to the resistivity. For the case of a potential V (r) = Vq 
at r < Ro and V(r) = at r > Ro, the estimation for this type of impurity contribution to the 
resistivity is[^] p ~ (h /4e 2 )«def^0' giving a negligible contribution to the resistivity when the radius 
of the potential Rq is of the order of interatomic distances and at a low concentration as above. 

3. Pentagonal and heptagonal defects in graphene lattice 

Exploring the 2D character and flexibility of this material, our idea is to propose a system 
which one or more sectors are excised from a graphene and the remainder is joined seamlessly. 
Removing a wedge from the graphene and identifying the edges a cone results j^]. In fact, the 
missed link of each carbon atom resting at the two edges of the remaining graphene sheet can 
be, in principle, covalently bounded. Particularly, considering the symmetry of a graphitic sheet 
and the Euler theorem, it can be shown that only five types of cones (incorporating one to five 



pentagons) can be made of a continuous graphene sheet [11, 12]. The motion of the charge carriers 
in an ideal conical graphene is equivalent to that of a massless Dirac particle in a gravitational field 
of a static particle of mass M in a (2 + 1 )D space-time. 

In the case of a cone with rid > 0, (the value rid, i n d = 1 , • • • , 5) is related to the conical angle 
Y , sin I = 1 — 2£ ) the deficit-angle induced by the conical singularity is given by 2^(1 — rid/6). The 
pentagonal defect can be presented as a pseudo-magnetic vortex at the apex of a graphitic cone, 
being the flux of the vortex related to the deficit angle of the cone (see Ref. []l3|]). Cones with a 
heptagon have negative curvature and are obtained by a insertion of a angular sector in the carbon 
sheet. Then if rid < 0, —rid counts the number of such sectors inserted into the graphene sheet. 
Our aim is, therefore, to see the influences that such a special graphene structure could induce on 
quasiparticle wavefunctions (spinors); surely, these influences may create new perspectives in the 
electronic transport properties, which are determined by the quasiparticles constrained to move on 
the conical surface. 
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To study the scattering of the carriers in graphene by topological defects we employ the anal- 
ogy between defects in condensed matter physics and in (2 + 1) -dimensional gravity Jl4| ] as far as 
possible. For example, the dynamics of charge carriers in an ideal conical graphene is equivalent 
to that of a massless Dirac particle in a gravitational field of a static point-like mass in a (2 + 1 )D 
space-time[15, [T^]. Specifically, we shall consider one of the simplest curved manifold, which is 
associated to the Schwarzschild solution in (2+ 1) dimensions: a space-time locally flat with global 
nontrivial properties. To describe this space-time we may use embedded coordinates r and 6 in the 
three-dimensional Euclidian space which extend over the complete range, O<r<oo 5 O<0< 2n, 
and describe a cone with the constraint z = \/ (cc~ 2 — 1)(* 2 +y 2 ), being the line element given by 
ds 2 = dt — a~ 2 dr — r 2 dQ 2 where a = 1 — 4GM Jl6[]. The attributes of the source are coded in the 
global properties of the locally flat variables. All the information lies in the non-trivial boundary 
conditions, which is important for the quantum scattering of graphene charge carriers by defects, 
like pentagons and heptagons. 

In the case of a topological defect in graphene, it is useful to change the gravitational term 4GM 
by the symbol /3, so that 271/3 (for < /3 < 1) gives the deficit of angle measuring the magnitude of 
the removed sector whereas —271(5 (for — oo < j8 < 0) accounts for the angle in excess associated to 
the insertion of a sector. The parameter /3 takes only discrete values because of the lattice symmetry 
of the graphene. 

Before analyzing the quantum mechanical scattering by conical defects (pentagons and hep- 
tagons) let us make a digression concerning scattering of the charge carriers in graphene as classical 
relativistic particles. The classical equation of motion, determined by relativistic geodesic equa- 
tion for the particles in a cone reads # + T^ [ „x a xP = 0, where the overdot indicates differentiation 
with respect to any convenient affine variable T that parametrizes the path x^(t) [[[(]]. The angle 
of scattering ±G) for the motion of the particles in a cone can be obtained by integration of the 
classical equations of motion and is given by Jl6| ] ± (0 = ±7i(a _1 — 1) = ±7t yAj , where ± refers 
to the side the charge carriers (current) trajectory pass around the defect. Note that the result above 
is valid for all values of j8 despite the sector was removed or inserted. The scattering angle above, 
presented in the embedded coordinate system measures the deflection of the asymptotic motion 
on cone projected onto x — y plane of the embedding three dimensional space. The result above 
suggests that a pentagon or heptagon may be used for deviating the planar current in graphene. 

To obtain the correct current deviations in graphene we have to solve the Dirac equation ( |EI| ) 
defined in a cone, say [|6p: ihy^Ea^D^ = 0, where D^ = d^ + \(0^ a bO ab , o ah = ^[y",/], and 
E a ^ is the dreibein in coordinates (t,r,8). The spin connection (0^- a b = —co^-ta mav t> e written in 
three dimensions as (0^ a b = £ a bcG>n c with £ a bc the Levi-Civita symbol as before [jig). The rotational 
invariance of the problem enables us to choose positive energy solutions that are simultaneously 
angular momentum eigenfunctions, with eigenvalue (n + j)Ti: 



u n (r)e- iEt l h = e ^-^ d [ "j^ ] e~ m l % , (3.1) 

<(r) 



where n = 0,±1,±2, ... . The solutions for E > are u^{r) = (£ n ) n J v (Kr), and u%(r) = 
(e„) n+1 J v (Kr). Here, J v is the Bessel function of order v = £ n /a(n + (1 =Fa)/2), n = 0, ±1, ±2 , ... , 
K = E/hvp(X,E > 0,£„ = ±1 and the same sign has to be chosen for the upper and lower compo- 
nents of the spinor u n {r). For 0<a<lor0<j8<l (remember that a = 1 — j8) we must choose 
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e n = sign (72 + (1 — CC)/2) = sign« (signO = 1) to have both components regular at the origin. The 
asymptotic form of the Bessel functions determines the phase shifts (they are identical for the upper 
and lower components) and are given by [|(]] : 

Sn = -^((l-a)n+(l-a)/2) = -^ r ^(n-i), (3.2) 

e n = sign(?i + (1 - a) /2) = sign(« + j8/2) . (3.3) 

The phase-shifts depend only on the number of sectors removed or inserted in the graphene sheet, 
accounted by a = 1 — j8. If — °° < j3 < 0, we need to be careful because e n = ±1 depending upon 
the value of (n + j3/2) (but the phase-shifts remains as above) and the phase-shifts depends only on 
the number of sectors (heptagons) inserted in the flat graphene sheet. In the presence of heptagons 
the carriers dynamics is identical to the that movement of the electrons in the gravitational field of 
a negative mass (although not possible in gravitation, this is feasible in the present context). 

Note that the phase shift ( |3~2| ) measures the deflection of the asymptotic motion on the cone 
projected onto x — y plane being qualitatively identical to the classical case discussed before. When 
there is a pentagon into the lattice and the fermionic current is constrained to pass around and 
sufficiently close to it such a current is scattered by the defect with an angle which depends only 
on the number of sectors removed in the graphene and on the side current passed. After passing by 
the pentagon the scattered current trajectories cross and yields an interference pattern. In the case 
of a heptagon, such a current is scattered but the trajectories diverge each other. 

4. Conclusions 

We have studied the scattering of graphene quasiparticles by topological defects like holes, 
pentagons and heptagons. We obtain the phase shift of the wave-function in all cases. For the case 
of holes, the main contribution concerns the s scattering and even in this case they do not change 
the resistivity of the sample, at least at low concentrations (like occurs to short range potential 
impurities). We realize that when the fermionic current is constrained to move near and around of 
pentagons and heptagons introduced in the lattice, it is scattered with an angle that depends only 
on the number of defects and on which side the current taken. Such a deviation may be determined 
by means of a Young-type experiment, through the interference pattern between the two currents 
scattered by the pentagon. In the case of a heptagon such a current is also scattered but it diverges 
from the defect. In addition, graphene would provide an appealing way to experimentally explore 
general relativity in two spatial dimensions since such effects are predicted by this theory [[l^, |l7| ]. 
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